"Fluid-solid phase-separation in hard-sphere mix- 
tures is unrelated to bond-percolation." 

In a recent letter, A. Buhot jjj proposes that entropy 
driven phase-separation in hard-core binary mixtures is 
directly related to a bond-percolation transition. In par- 
ticular, Buhot suggests that a phase- instability occurs 
when the coordination number rib, denned as: 

rib = Pi / gu(r)dr, (1) 

J CT ! <r<er ! (l+.R) 

is equal to zp c , where z is the coordination number of a 
particular crystal lattice, and p c is its bond-percolation 
threshold. Here pi is the number density of the larger 
particles, gu{r) is the radial distribution function of the 
larger particles, and R = a s /&i < 1 is the ratio of 
the diameters <7j. However, for binary hard-sphere mix- 
tures, calculations based on an accurate approximation 
to gu{r) demonstrate that n& varies widely along the 
phase-boundaries calculated directly by simulations, im- 
plying that bond-percolation is unrelated to the phase- 
separation in these systems. 

For highly asymmetric binary hard-sphere systems, Di- 
jkstra et. al. |2| conclusively demonstrated that an ef- 
fective one-component description based on a depletion 
potential picture quantitatively describes the fluid-solid 
transition. This in turn implies that the one-component 
description should give a fair representation of the ra- 
dial distribution function gu(r). Recent simulations |j] 
of the Asakura-Oosawa (AO) depletion potential |I| show 
that the Percus-Yevick (PY) approximation quantita- 
tively describes the pair-correlations along the fluid-solid 
transition line. In the inset of Fig. [I], the simple form: 
g u (r) = exp(-(3V de p(r)), where V dep (r) is the effective 
depletion potential, is compared to more the accurate 
PY integral equation results. Typically for packing frac- 
tion rji = irpicrf /6 < 0.25 along the phase boundaries this 
form gives near quantitative agreement for r < ai(l + R), 
which is not surprising since for small rji the potential is 
typically at least 2.5fceT along the fluid-solid transition 
line while the hard-core induced correlations are small so 
that the exponential form dominates. In fact, Buhot 's 
treatment of binary hard-spheres reduces exactly to this 
simple form but with an AO depletion potential which is 
valid only when 77; — ^ 0. If one replaces rj s with the rf s of 
small spheres in a reservoir kept at constant chemical po- 
tential, the correct form of the AO potential is recovered 
PL In Fig. [l], the metastable fluid-fluid and the stable 
fluid-solid phase lines taken from simulations Q are com- 
pared to lines of constant coordination number, which are 
calculated with eqn. ([!]) and gu{r) — exp (—/3Vd ep (r)), 
together with the depletion potential used for the simu- 
lations as well the simpler AO potential. The difference 
between the two is very small, implying that the coor- 
dination number is not very sensitive to the exact form 
of the depletion potential. Also included is the proposed 


bond-percolation induced fluid-solid line at rib = 2.4 de- 
rived from the approximation to gu(r) used by Buhot. 
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FIG. 1. Solid lines: the fluid-solid and fluid-fluid 
phase lines from the simulations of Dijkstra et. al. 
M for R = 1/30. Dotted lines and dash-dotted 
lines: the coordination number calculated as in the 
text using the AO potential and the potential from 
the simulations respectively. Long-dashed line: the 
n;, = 2.4 line proposed by Buhot 111. Inset: gu(r) for 
R = 0.2,7); = 0.25, rf s = 0.25 with AO potential. Solid 
line: from gu(r) — exp(— f3Vde P (r))\ dashed line: from 
PY integral equation. 

Clearly: (a) As expected, the approximation used by 
Buhot for gu(r) breaks down as rji increases, (b) The 
lines of constant coordination number are not related to 
either the fluid-fluid or the fluid-solid phase lines, im- 
plying that there is no direct relation between bond- 
percolation and phase-separation. The same results were 
found for other size ratios, and it is hard to see how more 
accurate approximations for gu (r) could change this pic- 
ture. 

The breakdown of the bond-percolation picture for this 
archetypical hard-core mixture model implies a similar 
breakdown for other, more complex, mixtures. The good 
agreement found at a few state points for hard-square 
systems probably results from either their 2-d nature, the 
imposed parallel symmetry, or the rather unusual pur- 
ported 2nd order fluid-solid transition. It does not imply 
that bond-percolation is generally relevant for fluid-solid 
phase-separation in binary hard-core mixtures. 
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